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Abstract 

Wood, due to its biological origin, has the capacity to interact with water. Sorption/desorption 
of moisture is accompanied with swelling/shrinkage and softening/hardening of its stiff¬ 
ness. The correct prediction of the behavior of wood components undergoing environmental 
loading requires that the moisture behavior and mechanical behavior of wood are consid¬ 
ered in a coupled manner. We propose a comprehensive framework using a fully coupled 
poromechanical approach, where its multiscale implementation provides the capacity to 
take into account, directly, the exact geometry of the wood cellular structure, using compu¬ 
tational homogenization. A hierarchical model is used to take into account the subcellular 
composite-like organization of the material. Such advanced modeling requires high reso¬ 
lution experimental data for the appropriate determination of inputs and for its validation. 
High-resolution x-ray tomography, digital image correlation, and neutron imaging are pre¬ 
sented as valuable methods to provide the required information. 
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1 Introduction 


Wood, an orthotropic cellular material, has the thought-provoking property of ad¬ 
sorbing water moleeules from its surrounding into its hierarehieal material struc¬ 
ture. Although known and employed for millenniums, wood in use and exposed to 
environmental loading is subjeet to varying moisture eontents (MCs) and eonse- 
quently exposed to a rieh pallet of moisture-indueed processes, possibly leading to 
wood degradation. 

Wood, eompared to other materials, is very hygroseopic. As water molecules attach 
themselves to the hydrophilie matrix in the cell walls, the induced fluid- solid-gas 
interaction forces result in a swelling of the eell walls. Due to cell geometrieal irreg¬ 
ularities and latewood and earlywood swelling incompatibilities, moisture-induced 
internal stresses originate and highly influenee the hygromechanieal behavior of 
wood as observed at the maeroseale. Wood is known as a material showing a high 
influenee of meehanical forees on water sorption. The interaction of the moisture 
and meehanical behavior of wood is best observed in swelling. Absorption of mois¬ 
ture in wood, in the hygroseopie range, that is, until around 30% MC mass per mass, 
results in swelling and reduced stiffness. The microseopie origin of this behavior 
lies at the seale of the eell wall. The eell wall material is eomposed in almost almost 
equal quantity of stiff eellulose microfibrils and a soft polymeric matrix. The hy¬ 
drophobic crystalline cellulose is surrounded by hydrophilie amorphous eellulose, 
immersed in a hydrophilic amorphous matrix, of hemicelluloses bound by lignin. 
The sorption of water moleeules in between the hydrophilie moleeules pushes the 
eonstituents apart, resulting in swelling and a reduction of stiffness of the matrix. 
The thin internal and external cell wall layers (i.e. S3 and Si) act as corsets due to 
the winding of the eellulose fibrils around the eell. In the eentral and, by far, thiek- 
est eell wall layer, namely S2, the eellulose micro fibrils are almost parallel to the 
longitudinal axis of the cells, although the presenee of an angle (called ’microfibril 
angle’ (MFA)) results in a helieoidal organization of the fibrils. Although the effect 
of varying the MFA from Oo to 30° on the transverse stiffness and swelling proper¬ 
ties of softwoods is small, during the sorption of moisture, the general orientation 
of the microfibrils in the S2 layer results in notable swelling in the transverse diree- 
tions of the cell and almost none along the longitudinal direetion. Furthermore, the 
helieoidal organization of the eellulose fibrils in the eell wall causes swelling to be 
more pronouneed normal to the eell wall than along the cell wall direetion. 

Although swelling originates from the eell wall level, its anisotropic nature seems 
to find its origin mainly at the eellular arehitecture level (see Figure 1 ). On a meso- 
seopie seale, wood consists mostly of longitudinal tracheid cells and also of radi¬ 
ally oriented ray eells, the latter representing 5 % of the wood volume. Across the 
growth ring, the thin-walled earlywood cells with large internal cavities, ealled ’lu¬ 
mens’, often gradually ehange to thick-walled latewood eells with small-sized lu¬ 
mens. The variation of the cellular strueture across the growth ring has been found 
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Figure 1. Tomographic image of spruce highlighting different anatomical features of soft¬ 
wood. 

responsible for the anisotropy of swelling and stiffness (Boutelje, 1962; Watanabe 
et al., 2000). Furthermore, recent findings on the isotropic swelling of latewood and 
anisotropic swelling of earlywood in the transverse plane indicate the possibility of 
latewood/earlywood interaction (Derome et al., 2011). 

At present, we do not completely understand the origins and play of water in the 
moisture, mechanical, and time-dependent processes in wood, not to mention their 
hysteresis (e.g. as investigated for paper by Derluyn et al., 2007) or its impact of 
durability (e.g. Viitanen et al., 2010). We present our current work in the under¬ 
standing of the role of water in the behavior of wood. A combined modeling/ ex¬ 
perimental approach is used to study further the role of water in these processes. A 
multiscale approach, using modeling and advanced experimental techniques, looks 
at the moisture-related properties of wood at different scales. The poroelastic mod¬ 
eling is used and upscaling procedures are developed from cell wall constituents to 
the scales of cell wall, cell, latewood and earlywood, growth ring, and wood tim¬ 
ber. To document and validate the modeling efforts, a comprehensive experimental 
plan has been established where synchroton X-ray tomographic datasets document 
the complex swelling/shrinkage behavior of wood cells, swelling at the growth ring 
scale is documented with digital image correlation (DIG), and finally, neutron pro¬ 
jections document changes in MC at the growth ring level. 


2 Modeling the water and mechanical behavior of wood 


We propose a comprehensive framework using a fully coupled poromechanical 
approach, where its multiscale implementation provides the capacity to take into 
account, directly, the exact geometry of the wood cellular structure, using compu- 
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tational homogenization. An hierarchical model is used to take into account the 
subcellular composite-like organization of the material. 


2 .1 Poromechanical approach 


A rigorous way of taking into account the interaction of fluids with the solid matrix 
in porous materials is based on the theory of poromechanics, introduced by Biot 
(1941). Within the context of thermodynamics of open porous continua, Coussy 
(2010) presented a general framework to formulate adequate constitutive equations 
for poroelastic behavior. The poromechanical approach starts from an energy view¬ 
point where stress/strain and chemical potential/MC are seen as pairs determining 
the energy content. With one equation expressing the energy content of the whole 
system, any change in stress field or chemical potential leads to a coupled response 
in both the strain and the MC. Thus, the poromechanical approach is suited to 
cover coupled effects between the mechanical and moisture response as observed 
in wood. A full description of the model described below is found in Carmeliet et 
al. (2013). 
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Figure 2. Poromechanical model results (continuous lines) fitted with experimental modu¬ 
lus of elasticity (black square) reported by Neuhaus (1981). 


In this coupled poromechanical approach, the free energy content, $ (J/m^), of a 
porous solid is described by the contributions of strain and moisture and the inter¬ 
action between strain and moisture. We can describe this system in terms of the 
applied fields, that is, stress and chemical potential, as follows 

= Sijdaij -|- ndp (1) 
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where eij (m/m) and aij (N/m^) denote the strain and stress tensors, respeetively; n 
(mol/m^) and m (J/mol) are the mole number density of the material and ehemieal 
potential of the adsorbed moisture, respeetively. 

The partial differentials of equation (1) yield two eonstitutive equations that de- 
seribe the eoupled meehanieal and moisture behavior of wood. The first of these 
two eonstitutive equations (equation ( 2 )) deseribes the strain (eij) with respeet to 
the stress (aij) and the ehemieal potential (/r). The seeond eonstitutive equation 
(equation (3)) gives the MC u with respeet to stress and ehemieal potential as fol¬ 
lows 


dsi — CijdcTj -\- Bidfi 

(2) 

du = -(Bidai + Mdp) 

(3) 


P 


where is the eomplianee (m^/Nm), Bi is the eoupling eoeffieient (mol/J), M is 
the moisture eapaeity(moP/m^/J), while p (kg/m^) and p (kg/mol) stand for volume 
density of the dried wood and molar mass of water, respeetively. 

Literature data were used to determine the introdueed parameters. Table 1 lists the 
fitted parameter. Figure 2 shows the eapaeity of the model to eapture the full range 
of eoupled moisture/meehanieal behavior of wood. 

Table 1 

Parameters fitted to the modef using data from Neuhaus (i981) for spruee. The seating 
vafue = f/12655 MPa is the eompfianee in the fongitudinaf direetion at zero moisture 
eontent_ 


ij 

^^3 — ^[0] 
^11 

3ij = ^ 

J[0] 

11 (LL) 

1 

3 

0.33 

22 (RR) 

13 

15 

0.14 

33 (TT) 

25 

16 

0.31 

12 (LR) 

-0.5 

17 

- 

13 (LT) 

-0.5 

29 

- 

23 (RT) 

-5.6 

38 

- 

44 (RT) 

222 

24 

0 

55 (LT) 

14 

14 

0 

66 (LR) 

17 

12 

0 
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2.2 Multiscale approach 


The interaction of the moisture and mechanical fields within a multiscale modeling 
framework requires a comprehensive treatment. The appropriate coupling between 
the scales is not straightforward. A rigorous multiscale framework in the context of 
a homogenization theory for the hygromechanical analysis of hierarchical cellular 
materials is thus sought. Given that the growth ring behavior is easily accessible 
and is considered here as the macroscale for wood, the effort is to develop an up- 
scaling technique from the cellular structure to the growth ring level. A detailed 
description of this work is found in Rafsanjani et al. (2012). We consider a problem 

(a) (b) 


r 




Figure 3. (a) Two-scale upscaling scheme and underlying REVs corresponding to different 
elements and (b) geometry of REVs. (REVs: representative elementary volumes.) 

of infinitesimally small deformations with moisture induced swelling for growth 
rings in softwood as a cellular porous solid within the framework of the compu¬ 
tational homogenization. Since the longitudinal dimensions of the wood cells are 
very long in comparison to the radial and tangential directions, the problem can be 
reduced to the analysis of the cross-section of the material, which is in a state of 
generalized plane strain. Verifying for separation of scales, the characteristic size of 
wood cells (microscale) is much smaller (fti40 /im) than the growth ring (macroscal 
) length scale (~3 mm). In the regions far away from the center of the stem of the 
tree, growth rings are periodically arranged in radial direction, which justifies an 
assumption on global periodicity of the timber. At the microscale, we assume local 
periodicity, where although wood cells have different morphologies corresponding 
to different relative ring positions, each cell repeats itself in its vicinity. Honeycomb 
representative elementary volumes (REVs) are selected for this current analysis. 
The concept of local and global periodicities in softwood, the proposed two-scale 
upscaling scheme, and the geometry of the REVs are illustrated in Eigure 3. 

We investigate the distribution of transverse anisotropy in swelling and mechanical 
behavior of wood from the cellular scale to the growth ring level by means of a finite 
element-based multiscale approach. The mechanical fields at the macroscopic level 
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Figure 4. (a) Lumen diameter and (b) cell wall thickness in radial and tangential directions 
as a function of the relative ring position, from SEM measurement on a spruce growth ring. 
(SEM: scanning electron microscopy.) 

(growth ring) are resolved through the ineorporation of the mierostruetural (eellu- 
lar strueture) response by the eomputational homogenization of different REVs of 
wood eells (honeyeomb REVs) seleeted from a morphologieal analysis of wood at 
the eellular seale (Derome et ah, 2012). Example of these measurements is shown 
in Eigure 4. 

These measurements are fitted to functions, which are used to predict the mechan¬ 
ical properties of unit cells along the growth ring using a master node homoge¬ 
nization technique. Eigures 5 and 6 compare the calculated results to experimental 
values. Eurther comparison of the results of the predicted swelling and elastic mod¬ 
uli at the growth ring level with experimental data, although not reported here, is 
also very satisfactory (Rafsanjani et ah, 2012). 



Eigure 5. Homogenized tangential and radial swelling coefficients for honeycomb unit ceiis 
(EEM finite element modeling, EXP experimental data from Derome et al. 2011). 


This multiscale approach achieves a more realistic characterization of the anisotropic 
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Figure 6. Homogenized tangential and radial elastic moduli for honeycomb unit cells (FEM 
finite element modeling, EXP experimental data from Earuggia and Perre, 2000). 

swelling and mechanical behavior of wood and, in particular, at describing the late- 
wood and earlywood interaction caused by the material heterogeneity at the cellular 
scale. 


2.3 Hierarchic micromechanical approach 


As mentioned above, wood interaction with moisture, although influence by all 
scales, stems from the moisture sorptive capacity of its constituents, that is, hemi- 
celluloses, celluloses, and lignin. Therefore, cell wall mechanical and moisture 
properties are not constant but depend on MC and the local spatial organization of 
the constituents. Even though the knowledge on the mechanical and moisture prop¬ 
erties of constituents is rather sparse, modeling can be used to validate fundamen¬ 
tal assumptions on constituents and their hierarchic organization. In a nutshell, in 
the wood cell wall, crystalline celluloses are surrounded by hemicelluloses, which 
are embedded in lignin, the latter binding together the cellulose/hemicellulose mi¬ 
crofibrils (Figure 7(a)). Microfibrils can be oriented unidirectionally or randomly in 



Eigure 7. (a) Constituents of a cell wall material and (b) cell wall of tracheid cells. ME: 
middle lamella; P: primary cell wall; CME: compound middle lamella. 

different cell wall layers. The angle between the longitudinal axis of a tracheid and 
the microfibrils is called microfibril angle (MFA). Five layers, namely, the middle 
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lamella (ML), primary cell wall (P), and the secondary ones (Si, S 2 , and S 3 ), build 
up the cell wall of the tracheid (Figure 7(b)). We focus on the hygromechanical 
behavior of cell walls. 



Figure 8. Calculated moisture- and microfibril-dependent elasticity of the cell wall com¬ 
pared with the results proposed by Persson (2000) at a moisture content of 12% (gray lines 
with circle markers). MFA: microfibril angle. 



Figure 9. Calculated swelling/shrinkage coefficients of the cell wall depending on the MFA 
and the thickness of the S 2 layer compared with the results proposed by Persson (2000) 
(gray lines with circle markers). MFA: microfibril angle. 


The layered nature of the cell walls, as well as the high aspect ratio of microfibrils, 
calls for a modeling approach analogous to the one used to model fiberreinforced 
composites. The first of the two main ingredients is the mixing rule, meaning the 
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derivation of effective layer properties from the constituent properties and, as such, 
a homogenization step. For this purpose, we employ a modified Halpin Tsai ap¬ 
proach with hygroexpansion. The second ingredient is the reduction of the five 
orthotropic layers to an equivalent single layer using classical laminate theory with 
the elastic properties of the constituents and the coefficients of hygroexpansion 
taken from Cave (1978) and Persson (2000). While our laminate approach follows 
Tsai and Hahn (1980), the mixture rules need to be discussed in detail. 


In comparative studies of mixture rules (Chamis, 1984), the Halpin-Tsai equations 
show the best predictive behavior. This is due to the fact that properties are not 
derived from first principles but are calculated using an empirical geometry factor 
for the spatial organization. This factor is estimated by fits to experimental data or 
detailed numerical models. For unidirectional (i.e. layers S 2 , and S 3 ) and ran¬ 
domly oriented (i.e. layers ML and P) materials, these factors are known (Halpin 
and Kardos, 1976). One peculiarity of micro fibrils is that the homogenization needs 
to be done in two steps, by combining first celluloses and hemicelluloses and then 
embedding this combination in lignin (Figure 7(a)). The two consecutive homog¬ 
enization steps for the coefficients of hygroexpansion ( a) were performed for the 
transversal isotropic Si layers as follows 


V}E) + 


(4) 


and for the longitudinal cell direction and for the directions perpendicular to the 
cells as follows 


+ (1 



( i + VfKA \ 

\V}Ej + V^El^j 


(5) 


where subscripts / and m denote the fiber phase and matrix phase, respectively, and 
V is the corresponding volume fraction. Due to their random-oriented character, the 
ML and P layers were subsequently averaged considering Young’s moduli in plane 
(El, E 2 ) given by 

_ 1 ^ aiEicos(^) + atE2Cos(^) 

- 3 g ■ E1 + E2 ^ ^ 


Our moisture-dependent analytical model is validated against a computationally 
expensive numeric two-scale finite element model proposed by Persson (2000) for 
one specific MC as a function of the MFA. Note that we used composite sections of 
ABAQUS for this purpose for simplicity and that the hygroexpansion can be treated 
similarly to thermal expansion. Figure 8 provides the calculated elastic properties of 
the cell wall, which are despite of small deviations in good agreement with Persson 
(2000) in all material directions. 
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To estimate the free swelling of a piece of cell wall, a block of two adjacent cell 
walls is calculated. The swelling coefficients in the main directions of the cell wall 
are given in Figure 9. It is not surprising that the swelling in thickness direction 
dominates the one seen in the other orientations. Moreover, the dependency on 
MFA is more pronounced in the longitudinal direction and in the direction parallel 
to the cell wall. Compared to Persson (2000), all tendencies can be confirmed, even 
though the absolute values perpendicular to the cell wall differ. The latter can be 
explained by different mixing rules used at the microfibril level. For consistency, 
we employ Halpin-Tsai as described above, while Persson (2000) uses a numeric 
homogenization scheme from a so-called base cell with adjusted geometry. 

In addition to provide the cell wall properties, this approach could be further used 
to study early/latewood for different cell wall and geometry properties in order 
to identify the factors dominating the swelling anisotropy of wood seen at the 
macroscale. 


3 Imaging methods 


The models presented above require high-resolution experimental data for the ap¬ 
propriate determination of inputs and for their validation. High-resolution tomog¬ 
raphy, Die, and neutron imaging are presented as valuable methods to provide the 
required information. 


3.1 Swelling of wood probed by phase-contrast X-ray synchroton tomography 


We investigated the three-dimensional (3D), microscopic, dimensional changes of 
spruce (Picea abies (L. Karst)) wood samples due to controlled steps of the ambient 
relative humidity (RH; see Figure 10). This study was performed at the wood cellu¬ 
lar scale by high-resolution synchroton radiation phase-contrast X-ray tomographic 
microscopy (srPCXTM), at the Tomcat beamline of the SLS of the Paul Scherrer 
Institute in Switzerland. 

Due to the low X-ray attenuation coefficient of wood, the phase-contrast method is 
used in order to reconstruct the information related to the X-ray index of refraction. 
In fact, the difference in density at the boundary between air and wood is a source 
of phase contrast. The reconstruction uses information from the coherent X-ray 
propagation through the specimen and interference of the X-ray photons coming 
out of the specimen. Phase contrast allows a more accurate determination of the 
boundaries between the components. 

Tomographic images of 500 pm x 500 pm x 8 pm samples, one of pure late- 


11 




(b) 


, 1 

Environmental < 

chamber \ 

1 

relative humidity 
sensors 





1 

1 

( 1 

sealed 

cap 


air inlet 


d 

window 

•'-"V l_ 


air outlet 


D 


sample 
and holder 








rotational 

stage 









Figure 10. (a) View of the experimental setup in the beamline and (b) schematic view of 
the environmental chamber. 

wood and the second of pure early wood, were taken as seen in Figure 11. Samples 
were let to achieve moisture equilibrium at five adsorption (starting at 25%RH to 
85%RH) and four desorption (down to 10%RH) steps plus a final adsorption step to 
25%RH. Affine registration of the dataset led to the identification of the orthotropic 
swelling strains, as shown in Figure 12, where the closure of the hysteretic loops 
demonstrates the reversibility of the process, and the start of a secondary loop, 
with the 25%-10%-25%RH steps, provides further evidence of the hysteretic be¬ 
havior. More details in the data processing that led to Figures 11 and 12 are found 
in Derome et al. (2011). 

Figure 12 shows that, for spruce latewood, swelling and shrinkage are found to be 
larger, more hysteretic, and more homomorphic than for earlywood. Furthermore, 
while latewood undergoes similar strains in the transverse directions, the early- 
wood radial strains are less than a third of the tangential strains. The less homo¬ 
morphic and smaller swelling/shrinkage of earlywood in radial direction is found 
to be caused by the presence of rays. The investigation of swelling/shrinkage strains 
was combined with dynamic vapor sorption (DVS) analysis, in which the samples 
were subjected exactly to the same RH conditions in adsorption and desorption. 
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The results show hysteresis in both latewood and earlywood, with MC values in a 
range between 4% and 18% within the same RH range, that is, between 10% and 
85%. 


A reeent investigation aimed at doeumenting the shrinkage of wood from green to 
dried state. The measurements were performed at the University Gent Center for X- 
ray Tomography, similar to the one deseribed in Massehaele et al. (2007). A sample 
of wood freshly taken from a live tree was subjeeted to several deereasing steps in 
RH. Using appropriately segmented 3D dataset, a eell wall volumetrie indieator is 
obtained by ealeulating the amount of material from the binary images at eaeh RH. 



Figure 11. Three-dimensional tomographie reconstruction of specimens: (a) latewood and 
(b) earlywood. 

In green state, a large pereentage of the sample eell lumens are filled with liquid 
water. As wet wood dried, free water leaves the lumen and the swelling pereent¬ 
age is obtained by eomparison of a volumetrie indieator in relation to the smallest 
volume measures. Similar to the previous work, a DVS analysis was performed to 
study the MC from the green to the dried states. The behavior in the first desorption. 
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and a subsequent sorption/desorption loop, is shown in Figure 13. 

The series of tridimensional dataset during swelling and shrinkage is then used as 
geometrieal input and for the validation of the multiseale modeling work presented 
above. 


(a) (b) 



RH (%| RH (%) 


Figure 12. Swelling strain measurements on pure latewood in (a) tangential and (b) radial 
directions and on pure earlywood in (c) tangential and (d) radial directions using microto¬ 
mography dataset. RH: relative humidity. 


3.2 Intraring variation of cross-sectional swelling by DIC 


In order to eompare the dimensional ehanges as generated by X-ray tomography on 
segmented earlywood and latewood samples with swelling at the maeroseale, we 
also investigated the moisture-indueed deformation on the growth ring seale using 
a eommereial DIC system (VIC2D; Correlated Solutions, USA). The noninvasive 
system enables to study the surfaee deformation of samples by eross-eorrelation 
of gray-value distributions. Therefore, a random gray-value speekle pattern was 
applied on the sample surfaees. The samples (40(R) x 40(T) x 5(L) mm^) were then 
eonditioned in a elimatie ehamber until they reaehed equilibrium, and an image of 
the sample was aequired. The images for eaeh RH steps were then analyzed with 
the image eorrelation software and refereneed to the dry-state image. 

The full deformation strain fields of a sample that underwent a MC ehange from 
oven dry (0%) to 16.84% are displayed in Figure 14. One ean elearly see that the 
growth ring strueture of the sample (Figure 14(a)) is elearly refleeted in the ra- 


14 












Figure 13. (a) Swelling indicator versus RH for first desorption from green state and ad¬ 
sorption, and second desorption loop, for spruce, (b) MC versus RH for the same sample. 
RH: relative humidity; MC: moisture content. 



Figure 14. Distribution of cross-sectional strains at moisture content 16.84% of the sample 
seen in (a) orthographic image, (b) tangential direction, and (c) radial direction. Scale of 
legend is swelling strain (%). 

dial strain distribution (Figure 14(c)), whereas the tangential strain field appears 
to be more or less homogeneous. In the radial strain distribution, the regions of 
higher deformation coincide with the latewood regions. Line extraction of radial 
and tangential strains is given in Figure 15. The average strain in radial direction 
is approximately 2.1%, whereas the tangential average strain is 4.8%, thus more 
than double. Furthermore, the more or less constant tangential swelling strain over 
the area of the sample results in a low deviation between maximum and minimum 
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strains (2.5%). For the radial strains, however, pronouneed differences can be seen 
for early wood and latewood resulting in the high deviation of 5.6%. These marked 
differences coincide with the radial cell wall thickness and thus density profile 
along the growth ring where the thin-walled earlywood cells result in low density 
and the brick-shaped thick-walled latewood cells in a significantly higher density 
(e.g. Derome et ah, 2011; Lanvermann et ah, in press). Thus, we can see that at the 
macroscale, latewood undergoes similar strains in tangential and radial directions, 
while earlywood undergoes low radial swelling strains compared to what happens 
in the tangential direction, confirming the findings presented above on the isolated 
tissues documented with X-ray tomography. 



Figure 15. Line extractions from the strain fields in radial and tangential directions. The 
radial strain distribution is displayed with solid line and the tangential one with a dashed 
line. 


3.3 Neutron radiography for moisture transport analysis 


Liquid water uptake in three softwood species (spruce, fir, and pine) in the radial, 
tangential, and longitudinal directions, is investigated using neutron radiography 
at the wood growth ring scale. The high sensitivity of neutron to hydrogen atoms 
enables an accurate determination of the change in the MC of wood. The analysis 
of the spatial and temporal change of water content distribution shows that liquid 
water transport has different characteristics depending on the direction of uptake. 
Figure 16 shows our custom-made experimental setup inside the thermal neutron 
radiography facility. A reference radiograph of the specimens at the initial state 
(dry or moist) is obtained after the sample is placed in the sample holder. Then, 
the water container is remotely moved up until the base of the sample is in contact 
with the free water surface. From the time of contact with water, radiographies are 
acquired at 60-s intervals, over a period of 60-180 min depending on the uptake 
rate of the different specimens. 

As a control reference, the total mass of intruded water during the uptake test is 
measured by weighing the mass of the specimens during the uptake test. To validate 
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Figure 16. Schematic overview of the experimental setup inside the NEUTRA beamline: 
(1) neutron source, (2) collimator, (3) boronated polyethelene blocks for shielding, (4) sam¬ 
ple installed on the sample holder, (5) detector, (6) remote control elevator, and (7) water 
container. 


dry 


moist 


width position (mm) 

Figure 17. Distribution of moisture content in fir specimens after 60 min of water uptake in 
longitudinal with dry or moist initial conditions (legend is in kg/m^). 

the moisture uptake eurves obtained from neutron radiographs, we also repeated a 
series of water uptake tests on the same samples in the laboratory using a preei- 
sion balanee. The detailed proeessing of the images is deseribed in Sedighi Gilani 
et al. (2012). Our results, see examples in Figures 17 to 19, show that latewood 
eells play a more significant role in water uptake than earlywood cells and that ray 
tracheids also contribute to liquid transport. Latewood tracheids possess smaller 
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width position (mm) 

Figure 18. Distribution of moisture content in fir specimens after 43 min of water uptake in 
tangential direction with dry or moist initial conditions. 

cell lumens and a smaller ratio of aspirated bordered pits than early wood eells that 
make them the preferential pathways for transport along the longitudinal direetion. 
The proeess of liquid uptake is different in the radial and tangential direetions as 
the path of the liquid is more intrieate, involving also the rays and requiring more 
often traversing pits. In tangential direetion, water uptake is oeeurring first in the 
latewood with a subsequent radial redistribution toward the earlywood. In radial 
direction, the growth ring boundary decreases the liquid transport rate, an indica¬ 
tion that a significant portion of the rays are interrupted at that loeation. All liquid 
transport is aeeompanied by sorption in the eell walls, but this sorption is most visi¬ 
ble during uptake in the radial direction. The moisture uptake rate in initially moist 
speeimens is seen to be higher. 

Neutron images can be proeessed to lead average MC along the height of the sam¬ 
ples, shown in Figure 20(a), (e), and (e). Furthermore, neutron images ean be reg¬ 
istered to lead swelling coefficients as shown in Figure 20(b), (d), and (f). 

The analysis of neutron images reveals interesting features about the proeess of 
liquid uptake and quantifieation and visualization of absorbed water, with high res¬ 
olution (0.07 kg/m^ in this experiment). We see a dependeney of the liquid water 
transport on the hierarehical strueture of softwood, orthotropie uptake direetions, 
and initial MC of the specimens. Higher uptake of the initially moist speeimens 
eould be attributed to their already modified (swollen) geometry or the existence 
of thin film of water on the mieroscopic pore surfaces, which increase the rate of 
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Figure 19. Distribution of moisture content in radial direction after 145 min of water uptake 
with dry or moist initial conditions. 

water absorption. 


4 Conclusion 


The proposed multiscale approach combines efficient computational tools with sup¬ 
porting swelling validation experiments at different spatial scales, as the mecha¬ 
nism of swelling links the moisture and mechanical behavior. The predicted swelling 
coefficients are compared to experimental data at the scales of the cell wall, of 
earlywood and latewood, and at macroscale. The agreement is good, predicting 
general trends accurately and giving a physical understanding of the observed phe¬ 
nomena. In addition, the complex structure of a natural material, such as wood, 
introduces a superposition of different influences that makes it difficult to eval¬ 
uate. The modeling framework has the flexibility to include different nonlinear 
and moisture-dependent constitutive material behavior at the microlevel and can 
be used for studying the influence of morphological, mechanical, and moisture ef¬ 
fects on effective hygromechanical properties, here, in softwoods and eventually in 
other wood-based composites. 

This comprehensive approach, combining modeling and advanced imaging, is ex- 
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Figure 20. Variation of moisture content along the height of pine samples during uptake 
in (a) longitudinal, (c) tangential, (e) radial directions, and (b, d, and f) corresponding 
swelling. 


pected to lead to a basic understanding of the cellular and hierarchic physics of 
wood in response to environmental stimuli, which could be a stepping stone for the 
understanding of other cellular or biological materials and the development of new 
biomimic materials with desired tailored properties. 
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